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Abstract 

There have been the approximate analytic solution [47] and several approximate analytic forms 
[18, 44, 45] of the growth factor D g for the general dark energy models with the constant values 
of its equation of state oj^e after Heath found the exact integral form of the solution of D g for 
the Universe including the cosmological constant or the curvature term. Recently, we obtained 
the exact analytic solutions of the growth factor for both oj& e = —1 or — | [48] and the general 
dark energy models with the constant equation of state Wdo [50] independently. We compare the 
exact analytic solution of D g with the other well known approximate solutions. We also prove 
that the analytic solutions for u>^ e = —1 or — ^ in Ref. [48] are the specific solutions of the exact 
solutions of the growth factor for general Lo^e models in Ref. [50] even though they look quite 
different. We scrutinize the issue of using the well known parameterizations of the growth index 
and its parameter given in Ref. [46] to obtain the growth factor for general dark energy models. 
We also investigate the possible extensions of the exact solution of D g to the time- varying Wd e for 
the comparison with observations. 



1 Introduction 



The analysis of the luminosity distance as a function of redshift obtained from distant Type la 
supernovae discovered that the present Universe is expanding at an accelerating rate [1, 2, 3]. 
One of the most popular solutions to this conundrum is introducing the so called "dark energy" 
(DE) which is the dominant energy contribution to the present energy of the Universe with its 
equation of state (EOS), u>de < — 3 (for example, [4]). The combined observations of the large scale 
structure (LSS) and of the cosmic microwave background (CMB) power spectra have confirmed 
the cosmic concordance (i.e. a flat universe with the present energy density contrast of the matter 
fl® n ~ 0.3 and with that of the dark energy fl° dc ~ 0.7) [5, 6, 7, 8, 9, 10]. 

Due to our ignorance of the nature of the dark energy, it is practical to use the EOS of the DE, 
u>dc to characterize it [11]. Moreover, u>de is the quantity constrained by cosmological observations 
[12, 13]. Among the excess of models, the cosmological constant A and a quintessence field are the 
most commonly proposed candidates for dark energy [14] . The former is characterized by code = — 1 
and the latter is a dynamical scalar field leading to a time dependent EOS, u;<j e (a). Also models 
with the constant uj^c = constant 7^—1 are important because the effects of the time varying 
ojde(o) can be predicted by interpolating between models with constant code [15, 16, 17, 18, 19, 20]. 

The origin of the current accelerating Universe is still in dispute (see for example, [21]). There 
are two major theories for this. One is the dark energy and the other is the modified theory of 
gravities (MG). However, MG are also able to be characterized by the effective EOS which is 
used for specifying DE models [22, 23]. Unfortunately, observations only probe the cosmological 
evolution of ujdc in an indirect way and there might be some ambiguities to differentiate DE 
with a specific MG model [24]. However, in most cases, while the two models give the same 
cosmic background expansion history H(a), they predict different growth rates for cosmic LSS 
[23, 25, 26, 27, 28, 29]. 

Thus, it is important to probe the accurate background expansion history of the Universe in 
order to constrain the EOS of the dark energy (i.e. its energy density, pd c ) precisely [12, 13]. 
Furthermore, the evolution of the matter density perturbation 5 m also depends on code [16, 30, 
31, 32, 33]. The formation of the LSS depends on the sound speed of the DE too [30, 32, 34, 35, 
36, 37, 38]. However, in general DE models including the quintessence, the sound speed of DE is 
close or equal to that of light and the DE is not able to cluster on the scales of galaxy clusters 
and below [30, 32]. Consequently, the DE only affects the matter power spectrum on large scales 
(> lOOMpc) [16]. Usually, the LSS measurements probe scales lOOkpc ~ lOOMpc and thus we 
may not need to worry about the effect of the growth of perturbation of DE when interpreting 
the LSS survey data. 

In sub-horizon scales (k 3> aH), all the matter density perturbation modes 5 m (k,a) grow 
uniformly because the dark energy do not cluster (^dc^dc *C ^ m <5m) and only the pressureless 
dark matter contributes to the gravitational potential. Thus, the effect of the existence of DE 
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appears only through the Hubble parameter, H(a) and one can use the linear growth factor D(a), 
defined by S m (k, a) oc <5# (fc)D(a). From the growth factor, the growth index / (sometimes it is 
called as "growth rate") is defined as / = dl ^ a) = Sl m (a) 7 [40, 41, 42]. In a flat universe, the 
growth factor is obtained in the integral form for the cosmological constant A [43]. This solution is 
widely used with the approximate analytic form [44] . This solution is even extended to the general 
dark energy models Wdc 7^ — 1 [45] by using the well known growth index parameter 7 (sometimes 
it is called as "growth index") given in the literature [46]. It is also known that the approximate 
analytic solution of D(a) is obtained in the general dark energy models with the constant uj^c 
[47]. However, the growth index f(a) = f2 m (a) 7 obtained from the approximate solution of 7 
in Ref. [46] causes the discrepancy with the correct one in certain DE models [48, 49, 50, 51]. 
The approximate solution of 7 is insensitive to the time ( i.e. a) and this cause the problem in 
interpretation of other quantities obtained from this approximate 7. Also the known approximate 
analytic solution of D(a) in Ref. [47] does not show the proper behavior of D(a) in some values 
of and/or some DE models [51]. 

One may use the definition of f(a) with the assumption of slowly varying or constant 7 as in 
Ref. [46] to obtain the values of D{a) and f(a) close to the real D(a) and f(a) for certain values 
of LUde and [46]. However, D(a) and f(a) obtained from this approximate value of 7 show the 
discrepancies with the correct ones in some values of oo^c and/or fi^ [48, 50, 51]. Thus, if one 
does any analysis related to the growth factor and/or growth index with the assumption that 7 
is almost constant in general DE or MG models, then the results obtained from this 7 are not 
reliable [51]. We have currently available data for /(a) at various redshifts with the large degree 
of uncertainty though [9, 52, 53, 54, 55, 56, 57]. 

In what follows, we analyze in detail the recently obtained exact analytic solution of the growth 
factor D(a) with the general constant uj^c dark energy in a flat universe [48, 49, 50, 51]. We note 
that the well known analytic solution of D(a) in Ref. [47] is the approximate solution which shows 
the different behaviors of both D(a) and /(a) from the exact ones for some DE models. We do 
confirm that the exact analytic solutions of the growth factor with to^e = — 1 and — | obtained in 
Ref. [48] are the specific solutions of the exact solution of D{a) for general oj& e given in Ref. [50] 
even though they look quite different. In Sec. 3, we compare the cosmological evolution of D(a) 
obtained from the well known approximate analytic forms of it with those of the exact analytic 
solution D(a). We also compare the values of f(a) from these two solutions. We analyze the 
problem of the well known approximation of 7 given in Ref. [46] when one use it to obtain D(a) 
and f(a) in Sec. 4. We investigate D(a) and /(a) with a specific parametrization of code and its 
applications to observations in Sec. 5. We reach our conclusions in Sec. 6. 
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2 Sub-horizon scale growth factor 



We use the flat Friedmann- Robertson- Walker universe to probe the sub-horizon scale linear density 
perturbations of matter 5 m in the matter dominated epoch, 
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where u; dc is the equation of state (EOS) of dark energy, p cr is the critical energy density, p m 
and pdc are the energy densities of the matter and the dark energy, respectively. We consider 
the constant u>dc and set the present scale factor ao = 1. At sub- horizon scales (k » aH), all 
interesting modes of the matter density perturbation S m (k,a) grow uniformly as long as the dark 
energy do not cluster [16, 31, 32]. It means that we only consider the matter perturbation in 
Poisson equation in this scale. Thus, the growth factor D g (a) is defined as 



S m (k,a) oc 5 H {k)D g (a) , 



(2.3) 



where 5nik) is the scalar amplitude at the horizon crossing generated during the cosmic inflation. 
The alternative definition of the D g are also commonly used (see for example, [39]) 



S m (k,a) = 5 (k)D g (a) , 



(2.4) 



where So(k) is the present density contrast. Then we obtain the evolution equation of D(a) from 
the linear density perturbation equations [40, 58], 
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We use D in Eq. (2.5) instead of D g because the general solution of Eq. (2.5) does not guarantee 
that D is the growing mode solution. We are able to find the exact analytic solution of D(a) 
for any value of the constant oj& c [48, 49, 50]. After replacing new parameters Y = Qa 3u)dc and 

Q = 



in Eq. (2.5), we obtain 
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We replace a trial solution D(Y) = cY a B(Y) into Eq (2.6) to get 
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W m = 0.2 


n» m = 0.3 


W m = 0.4 


Cswl 


Csw2 


Cswl 


Csw2 


Cswl 


Csw2 


-1/3 


0.305329 


0.370645 


0.484579 


0.546607 


0.723051 


0.783376 


-0.8 


0.563274 


0.568201 


0.704043 


0.707640 


0.845712 


0.848492 


-1.0 


0.630418 


0.631792 


0.754267 


0.755227 


0.873819 


0.874533 


-1.2 


0.680498 


0.680868 


0.790355 


0.790606 


0.893532 


0.893715 



Table 1: c sw i and c SW 2 are the values of the coefficient c sw obtained from the two initial conditions 



of the growing mode solution, D(a,i) ~ a« and ~ 1, respectively. 



The above equation becomes the so called "hypergeometric" equation when we replace X = —Y, 
which has the complete solution [59] 
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where F is the hypergeometric function. Thus, the exact analytic solution of the above equation 
(2.5) is 
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(2.9) 



D(a) in Eq. (2.9) is just the general solution of the second order differential equation (2.5). Thus 
it does not have any physical meaning yet. It may represent the growing mode, the decaying 
mode or none of them before we choose the integral constants c\ and ci- If we want to have 
the correct growing mode solution from the above analytic solution, then this solution should 
follow the behavior of growing mode solution at an early epoch, say ai ~ 0.1. In other words, the 
coefficients of the general solution should be fixed by using the initial conditions of the growth 
factor, 
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After we fix the coefficients from the initial conditions, we are able to determine the growth factor 
D g (a) from the general form of solution D(a) in Eq. (2.9). If one want to obtain the decaying 
mode solution D^a) from Eq. (2.9), then one need to adopt the decaying mode initial conditions 
to obtain the correct coefficients 
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Now we compare the exact growth factor in Eq. (2.9) with the well known approximate growing 
mode solution [47], 
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We rewrite the second term in Eq. (2.9) using the linear transformation formula of hypergeometric 
function [59], 
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Thus, D™(a) in Eq. (2.12) becomes equal to D g (a) in Eq. (2.9) if and only if 
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(2.15) 



where we use the notations that c± g and C2 9 are the values of coefficients c\ and C2 obtained 
from the growing mode initial conditions in Eq. (2.10). D sw (a) given in Eq. (2.12) have several 
problems. First, D sw (a) contains only one integral constant c sw even though it is obtained from 
the second order differential equation. This problem might be solved if the two integral constants 
in the general solution (2.9) satisfy the conditions in Eqs. (2.14) and (2.15) simultaneously. In 



other word, D s 



would be an exact growing mode solution if the values of c sw obtained from 



both initial conditions (2.10) are the same. We denote that c sw i and c sw 2 are the values of c s 



obtained from the growing mode initial conditions D„ 



~ dj and 



~ 1, respectively. As 



we show in Tab. 1, c sw i and c sw 2 show discrepancies for the different models. As Ude decreases, 
the difference between the two coefficients also decreases. The same effects happen when is 
big. Thus, D sw is a good approximate solution for the small value of u>dc and the big value of f^. 



One may suspect that this discrepancy between c sw i and c sw 2 might be due to the choices of 
initial conditions. We investigate this as follows. The exact values of initial conditions can be 
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c swl 


Csw2 


Awi(i) 


A w2 (i) 




fsw{ a i) 


f(ai) 


-0.4 


-0.0486 


0.589752 


0.606695 


0.624667 


0.642613 


0.631619 


0.924788 


0.951356 


-0.8 


-0.0011 


0.705845 


0.708663 


0.733052 


0.735979 


0.734198 


0.994917 


0.998890 


-1.0 


0.0000 


0.754267 


0.755227 


0.779311 


0.780303 


0.779699 


0.998730 


1.000000 



Table 2: £, c sw i, c sw 2, D sw i(a = 1), and D sw 2(a = 1) for the different values of w^e when we 
choose = 0.3 and cij = 0.1. D sw i(l) and D SW 2(1) are the present values of the growth factors 
when we choose c sw i and c sw 2, respectively. D g {l) is the present value of the exact growth factor. 
/sw(oi) and f{cii) correspond to the initial values of the growth index obtained from D sw and D g , 
respectively. 



obtained numerically from Eq. (2.5), 

D g {a l ) = a\ + t and ^[ = (l+£)a|, (2.16) 

where £ indicates the deviation of the growth factor from the linear growth Oj at the initial epoch. 
We show the magnitudes of these £s for the different DE models in Tab. 2. In this table, we 
choose = 0.3 and cij = 0.1. As Wdc decreases, the value of £ also decreases because there 
is more matter component at the initial epoch a\ for the smaller values of uo& c . Thus, D g (oi) is 
close to ai for the smaller value of uo^c- We also shows the values of c sw i and c sw 2 for different 
DE models with initial conditions in Eq. (2.16). If we compare c sw i and c sw 2 values in Tab. 1 
with those in Tab. 2, then we find that the discrepancies in the two values are not removed even 
with the exact values of initial conditions. Thus, the deviations in c sw i and c sw 2 are the intrinsic 
problem of the solution D sw and irrelevant to the accuracies of initial conditions. We define D sw i 
and D sw 2 as the growth factor obtained from Eq. (2.12) when we choose the coefficient c sw as 
c sw i and c sw 2, respectively. The discrepancies in the present values of the growth factor D sw i and 
D SW 2 are decreased as oo^e is decreased. D sw (a) is a good approximate solution for the small values 
of LOdc and big values of fi^. We show this in Fig. 1. When lo^c is bigger than —1 and fi^ is 
small, -Dswi and D SW 2 show the big discrepancies with the correct D g as shown in the left panel of 
Fig. 1. The dashed, solid and dotted lines correspond to D sw 2, D g , and D sw ±, respectively when 
ujde = —0.4 and = 0.1. However, when ui^e is close and smaller than —1 and is big, both 
D sw i and D sw 2 are very close to D g as shown in the right panel of Fig. 1 where we use Wde = — 1-0 
and = 0.3. 



Second, one is not able to separate the growing mode solution and the decaying mode solution 
in D(a). So far, there are only two known possible cases for separating the two modes when 
u>de = — | and —1 [60]. This separation is impossible for general values of u>de as shown above. 
Thus, D sw is not the correct growing mode solution. D(a) itself in Eq. (2.9) is the solution of the 
equation (2.5), and this solution cannot be separated as the growing mode solution or the decaying 
one. As we explained in the above, after we obtain the general solution of the equation (2.5), the 
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Figure 1: Evolutions of the growth factors Ds. a) D sw 2, D g , and D sw i (from top to bottom) for 
ujde = —0.4 and = 0.1. b) D sw 2, D g , and A> w i (from top to bottom) for L0^ e = —1.0 and 
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0.3. 



solution D{a) can be interpreted as the growing or the decaying mode solution by applying the 
appropriate initial conditions given in Eqs. (2.10) and (2.11) to the general solution in Eq. (2.9). 

The third problem is related to the growth index / = ■ If we choose D sw as the growing 
mode solution, then the growth index becomes 
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Since D sw has only one coefficient, the growth index obtained from D sw is independent of c, 
we choose the exact values of initial conditions given in Eq. (2.16), then the value of the growth 
index at the initial epoch will become 

/(a = a i ) = l + £ + 0(£ 3 ) + --- • (2.18) 

Therefore, / sw (a = dj) is not same as /(a^) given in Eq. (2.18). This problem also happens when 
we choose the approximate initial conditions (2.10). Thus, the value of the growth index obtained 
from D sw shows the intrinsic discrepancies with that obtained from the correct growth factor D g 
as shown in Tab. 2. We find that the present value of / for = 0.4 should be close to 0.6 
independent of Wde and thus Fig. 3 in Ref. [47] is incorrect. 

Recently, we have also obtained the exact analytic solution of D{a) for u^e = — 1 [48]. There 
we have found that the solution of D g for uj^ = — 1 is given by 
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The form of Dg{a) looks quite different from D g (a) in Eq. (2. { 
general solution D{a) becomes 
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where we use the relation F 
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y/1 + Y in the second and the third 
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equalities [59]. Thus, the solution Dg(a) given in Eq. (2.19) is one of the particular solutions of 
D (a) when u^c = — 1- We are also able to obtain the particular solution of D(a) when ui^e = — j 
by using the same relation. 
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where arctanh is the inverse hyperbolic tangent function. 



3 Comparison with known approximate solutions 



There are several well known approximate analytic forms of the growth factor [18, 44, 45]. For 
the cosmological constant (i.e. uj^ c = — 1), the well known approximate form of the growth factor 
at present is given by [44] 

-i -i 
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where 0° is the present value of the energy density contrast {p%/Pc T ) °f the cosmological constant 
A. One is not able to obtain the growth index from the D® pt because it is a constant. Thus, the 
approximate analytic form of the growth index /i a hav is given separately in Ref. [42] : 
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The above solution is generalized to any value of a in Ref. [18] : 
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We compare this solution D cpt with the exact one in Eq. (2.20). 
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Figure 2: a) Cosmological evolutions of the growth factor D g and D cpt (from top to bottom) for 
= 0.3. b) Evolutions of the growth index / and f cp t (from top to bottom) for same f2^. 



We show the cosmological evolutions of D cpt and f cp t as compared to the exact solutions D g 
and / when = 0.3 in Fig. 2. The analytic approximate solution D cp t is very close to the 
correct growth factor D g as shown in the left panel of Fig. 2. The solid line describes the D g 
and the dashed line depicts the D cpt . The error in D cpt is less than 1 % for any a. / cp t(«) is also 
obtained from D cpt in Eq. (3.3). We compare the evolution of / cpt with that of / in the right 
panel of Fig. 2. The solid line describes the correct / obtained from D g . The dotted line shows 
the evolution / cpt obtained from D cpt . The error of / cp t at present is about 4 % only when we use 

= 0.3. 



There is also another approximate analytic solution D\, p for the general values of Wde [18, 
45]. This solution is obtained from the well known parametrization of the growth index and its 
parameter in Ref. [46], 

/ = ^ £ = ^m(a)^, (3.4) 
a ma 

3(l-a; de ) , 3 (1 - a; de )(l - jf*) 

where 7 WS ~ « — (1 — U m (a)j . (3.5) 

' 5 - 6wd e 125 (1 - 2^)3 v mv " v ' 

The approximate growth factor D^ p is known as the extension of D cpt in Eq. (3.3) for the general 
Wde an d is given by [45] 
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de 



(3.6) 



where 7^, is the approximate form of the growth index parameter by choosing f2 m (a) = in 
Eq. (3.5) and A is the fitting coefficient in Ref. [18, 45] 

^ 3(1-^) 3 (1-^(1 -3^) 
7ws ~ 5-6u dc + 125 (i_^)3 11 " mj ' V- l > 
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Figure 3: a) Cosmological evolutions of the growth factor D g and D^ p (solid and dashed lines, 
respectively) for lu Ac = —0.8 and = 0.3. b) Evolutions of / and /b p (with same notation) for 
same uj de and f2^. 



A ~ 1.742 + 3.343u; dc + 1.615u4 when to de > -1 
0.28 

~ 0.3 when w d e < — 1 . 

uj Ac + 0.08 



(3.8) 
(3.9) 



Again the above solution D® is generalized as [18] 



5Q m (a) 
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Even though the value of Db P at an Y epoch a is very close to that of D g , its evolution behavior 
is quite different from that of D g . We show this in Fig. 3. In the left panel of Fig. 3, the solid 
line corresponds to D g and the dashed line depicts D^ p with A given in Eq. (3.8) for uj de = —0.8 
and = 0.3. Even though the error in D^ p at present is about 3 %, there is the discrepancy in 
cosmological evolution behaviors of D g and Dbp- This discrepancy is clear when we compare the 
growth index / and /b p as shown in the right panel of Fig. 3. The solid line describes / and /b p 
is depicted by the dashed line. The error in /b p is more than 20 % at present. 

As we show in Figs. 2 and 3, one might be able to obtain the value of the growth factor with 
small error from the approximate analytic solutions of the growth factor. However, one needs to 
pay attention when one considers the growth index. Especially, might not be used to compare 
with observations because of the incorrect behavior of /b p obtained from the approximate analytic 
solution Z?b P for some DE models. 
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Figure 4: a) Cosmological evolutions of the growth factors D g and D ws (from top to bottom) for 
code = —0.4 and = 0.2. b) Evolutions of D g and D ws (from top to bottom) for cjde = —1-0 and 
fl° m = 0.3. 



4 Problems of D(a) and f(a) obtained from the approximate 7 

Nowadays, the parameterizations of the growth index / given in Eq. (3.4) and its parameter 7 in 
Eq. (3.5) are commonly used to derive the growth factor for general DE and/or MG [27, 46] 



D ws (a) ~ aexp 



f ^(tt m (a') 7ws -l 
a' V 



(4.1) 



Eq. (4.1) is slightly changed from the original D ws in Ref. [46] which is normalized to the growth 
factor at the present epoch [27]. It is known that the growth factor D ws and the growth index 
/ ws obtained from these parameterizations are very close to the correct values of D g and / [27]. 
However, this is true for small values of u>de and/or for big values of ft^. The evolution of the 
growth factor D ws based on this parameter ~f ws shows the discrepancy with the correct D g for the 
values of uj& e bigger than —1 [51]. Both growth factor and growth index are important to probe 
the power spectrum, the evolution of the mass function, the large scale structure statistics like in 
the gravitational lensing, the peculiar velocity, and so on. Thus, we need the accurate values of 
both D g and / to probe the DE and/or MG models correctly. 



We show the deviation of D ws from D g for the different values of cu^e in Fig. 4. In the left panel 
of Fig. 4, the dashed line describes the evolution of the growth factor D ws obtained from Eqs. 
(4.1) and (3.5) when Wd e = —0.4 and = 0.2. The solid line depicts the correct D evolution for 
the same values of Wde and Q^. The error of the current value of the growth factor is about 16%. 
However, the error in D ws is decreased for the smaller value of Wde and the bigger value of f2^. 
We show this in the right panel of Fig. 4. The solid and dashed lines correspond to D g and D ws , 
respectively in the right panel of Fig. 4 when we choose Wde = —1-0 and £1^ = 0.3. In this case, 
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0.2 0.4 0.6 0.8 1.0 0.90 0.92 0.94 0.96 0.98 1.00 

a a 

Figure 5: a) Cosmological evolutions of the growth index / and / ws (from top to bottom) for 
u>de = —0.4 and = 0.2. b) Evolutions of / ws and / (from top to bottom) for u>d e = — 1.0 and 
n° m = 0.3. 

both growth factors are almost identical for the entire of a. Thus, if one investigates the models 
based on the cosmological concordance model with slowly changing code, then the fitting formulae 
given in Eq. (4.1) and Eq. (3.5) are good ones. However, one should be careful to use these fitting 
formulae when one probes any model with time- varying uj^c and/or with different cosmological 
parameters from the concordance model. 

We also show the cosmological evolutions of / ws and / in Fig. 5. In the left panel of Fig. 5, the 
solid line depicts the growth index / when u>d e = —0.4 and $7^ = 0.2 and the dashed line shows 
/ ws for the same parameters. The deviation of / ws from the correct / is very small except at early 
epoch. We also describe / ws and / for the different values of Wdc and in the right panel of Fig. 
5. The dashed and the solid lines correspond to / ws and /, respectively when ujd e = —1.0 and 

= 0-3. Again we find that the error in / ws is decreased when u>de decreases and fi^ increases. 



Even though one is able to obtain the very close values of physical quantities like D and / to 
the correct ones from the fitting formulae, one should be careful for extending those formulae in 
general cases. Especially, one needs to pay more attention when one adopts them to differentiate 
DE with MG. We obtain the analytic solutions for D g and / that are exact for any DE model. 
And these solutions give the exact theoretical values of observable quantities. However, the exact 
solutions are limited to the constant values of uJde- Thus, we need to investigate the generalization 
of the solutions to more general cases including time- varying code- We will explain the possible 
extensions of them in the following section. 
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Figure 6: A specific model of Eq. (5.1). 



5 Applications for D g (a) and f(a) to time- varying u ( \ c 

Even though the growth factor obtained in Eq. (2.9) is only true for the constant u^ei we are able 
to apply this solution to the time-varying uj^e by interpolating between models with constant uj^lc 
[17, 18, 19, 20]. For this purpose, we choose the sum of the step functions 6 of u^ e (a) to probe 
the evolutions of D g (a) and f(a), 



where Wde(j) is the arbitrary value we need to fit from the background evolution observations. We 
use a specific model of this, to^ p = — 0.89(a) — 0.16 (a— 0.6) — 0.16 (a— 0.7), for the demonstration as 
shown in Fig. 6. The values of LOde(j) and a,j are related to the values of the uide parameterization 
which produce the proper background evolution like H(a) [19]. Also, one is able to extend this 
parametrization to more general cases by putting more steps and/or different values of u>d e (j). 



The advantages of this parametrization of u^c are the followings. Even though the EOS is a 
discontinuous function of a (i.e. z), the physical quantities like H(a), D g (a), and f(a) obtained 
from this wf ep are smooth functions [17]. We are able to obtain the smooth functions D g (a) and 
f(a) by solving for the proper values of c\ and C2 in Eq. (2.9) at each interval. This is shown 
in Tab. 3. Also the observations constrain the physical quantities in the specific interval of a. 
Thus, the parameterization of u)& e in Eq. (5.1) is a good one to probe the properties of u>d e when 
compared to the observations. 




(5.1) 
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0.1 < a < 0.6 


0.6 < a < 0.7 


0.7 < a < 1.0 


D g (0.1) 




D g {0.6) 


"atIo.6 


A,(0.7) 


"TErlo.7 


0.1 


1 


0.530531 


0.661798 


0.592623 


0.579974 


Cl 




Cl 




Cl 


C2 


1.09339 


-1.39169 


1.04914 


-1.07846 


1.01286 


-0.872337 



Table 3: ci and C2 are the values of the coefficients obtained from D{a 
each interval. 



and 



dD{a) 
da 



at 




(1.5 0.6 0.7 0.8 0.9 1.0 0.5 0.6 0.7 0.8 0.9 1.0 

a a 



Figure 7: a) Cosmological evolutions of D g (a) for uj& c = —1.0, c^ e ep , and —0.8 (from top to 
bottom) when = 0.3. b) Evolutions of /(a) for uj^ c = —1.0, —0.8, and ul ep (from top to 
bottom) for ^ = 0.3. 



We show the cosmological evolutions of D g and / in Fig. 7 for the different DE models. The 
evolutions of the growth factors D g for uj^c = —1-0, w^ e ep , and —0.8 are described as dotted, 
solid, and dashed lines, respectively in the left panel of Fig. 7. The present values of D g are 
(0.7797,0.7328,0.7327) for w do = -1.0, uj^ p ', and -0.8, respectively. The evolutions of the D g 
for uj^ p and Ud e = —0.8 are quite similar to each other because of the specific choices for values 
of lo^c P in Eq. (5.1). If we have the shorter period of aj=\ for Udc = —0.8 in then we 

may have the different evolution of D g from the different choices of u^ p . Also the cosmological 
evolutions of / are depicted as the dotted, solid, and dashed lines for lu^c = —1-0, —0.8, and wj e ep , 
respectively in the right panel of Fig. 7. We obtain the present values of /, (0.5128, 0.5084, 0.4974) 
when u>de = —1-0, —0.8, and wj e ep , respectively. Thus, we obtain very interesting features of D g 
and / from these DE models. Even though the present values of uj^e = —1-0 and Wj* ep are equal, 
the evolution behaviors of D g and / are quite different for these two models as shown in Fig. 7. 
D(a = 1) values are different by as large as 6 % and the difference in /(a = 1) is about 3 %. 
Thus, we may have a good chance to tell whether ui<ie is a constant or not by investigating D g {a) 
and f{a) at different a intervals. 
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6 Conclusion and discussion 



We have analyzed the properties of the exact analytic solution of sub-horizon scale matter density 
perturbation (i.e. growth factor) for the general dark energy models with its equation of state tu^c 
being constant. From the comparison of this solution D g with the well known approximate analytic 
solution D™, we have found that Dg W is a good approximate solution of D g for the concordance 
model. D g can be expressed with the slightly different functional forms for the specific values of 
code- Especially, we have explicitly shown the alternative forms of D g when u>& c = — 1 and — |. 
The two solutions in Refs. [48] and [50] are equivalent when oj^c = —1 or — | even though they 
look quite different. 

We have investigated the several well known approximate analytic forms of the growth factor. 
D cpt is the one with the dark energy being the cosmological constant and D^ p is the extension of 
D cpt for the general dark energy models with constant lo^c- /cpt and /bp are the growth indices 
obtained from D cpt and -Dbp, respectively. When the dark energy is the cosmological constant, 
D cp t and f cp t are very close to the correct D g and /. However, D^ p and /b p show the discrepancies 
with the correct D g and / for some dark energy models. Especially, the error in /b p for tu^e = —0.8 
and = 0.3 is as large as 20% at present. 

We have also scrutinized the validity of the well known parameterizations of the growth index 
/ ws and its parameter 7 WS . The growth factor D ws and the growth index 7 WS are very close to the 
correct ones for the cosmological concordance model. However, the error in the present value of 
the growth factor D ws obtained from these parameterizations is about 16% for oj^c = —0.4 and 

n° m = o.2. 

The approximate analytic solution D sw , the approximate analytic forms D cpt and -Dbp, and 
D ws obtained from the parameterization of 7 WS are good approximate solutions of the exact D g 
for the concordance model. However, all of them show some discrepancies with the correct D g 
for some DE models and/or 17^ values. Thus, one needs to be very careful when one extends the 
approximate solutions to the general models and/or other cosmological parameters. 

Even though we have obtained the exact analytic solution of D g for the general DE models, 
this solution is limited to the constant Wdc models. Thus, the applications of this solution to the 
real observations are very limited. However, we can apply this solution to the more general cases 
like the time- varying uj^c by interpolating between models with constant u>de- We have found that 
D g and / obtained from the constant uj^c and the time- varying one are quite different even though 
we have the same values of u;<ic s at present. If we are able to obtain a good constraint on lj^c from 
the cosmological background evolution observations, then we will be able to constrain D g and / 
very accurately. Thus, the exact analytic solution of D g can be used as the very useful tool for 
the interpretation of LSS survey data. 
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